ADP-03-123/T560 



Comparison of \Q\ = 1 and \Q\ — 2 Gauge-Field Configurations 
on the Lattice Four- Torus 

Sundance O. Bilson-ThompsonQ Derek B. LcinwcbcrQ and Anthony G. William^] 

CSSM Lattice Collaboration, Special Research Centre for the Subatomic 
Structure of Matter and The Department of Physics and Mathematical Physics, 
University of Adelaide, Adelaide SA 5005, Australia 

Gerald V. Dunnc| 

Department of Physics, University of Connecticut, Storrs CT 06269, USA 
(Dated: February 1, 2008) 

It is known that exactly self-dual gauge- field configurations with topological charge \Q\ = 1 cannot 
exist on the untwisted continuum 4-torus. We explore the manifestation of this remarkable fact on 
the lattice 4-torus for SU(3) using advanced techniques for controlling lattice discretization errors, 
extending earlier work of De Forcrand et al for SU(2). We identify three distinct signals for the 
instability of \Q\ = 1 configurations, and show that these signals manifest themselves early in the 
cooling process, long before the would-be instanton has shrunk to a size comparable to the lattice 
discretization threshold. These signals do not appear for the individual instantons which make up 
our |Q| = 2 configurations. This indicates that these signals reflect the truly global nature of the 
instability, rather than the local discretization effects which cause the eventual disappearance of the 
would-be single instanton. Monte-Carlo generated SU(3) gauge field configurations are cooled to 
the self-dual limit using an C2(a 4 )-improved gauge action chosen to have small but positive 0(a G ) 
errors. This choice prevents lattice discretization errors from destroying instantons provided their 
size exceeds the dislocation threshold of the cooling algorithm. Lattice discretization errors are 
evaluated by comparing the 0{ a 4 ) -improved gauge-field action with an 0(a 4 )-improved action con- 
structed from the square of an 0(a 4 )-improved lattice field-strength tensor, thus having different 
C(a 6 ) discretization errors. The number of action-density peaks, the instanton size and the topolog- 
ical charge of configurations is monitored. We observe a fluctuation in the total topological charge 
of \Q\ = 1 configurations, and demonstrate that the onset of this unusual behavior corresponds with 
the disappearance of multiple-peaks in the action density. At the same time discretization errors 
are minimal. 
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I. INTRODUCTION for multi-instantons are known on R 4 (and on S 4 ) via 

Oh, the ADHM construction 0, @ , no such concrete results 



are known on the continuum four-torus. The only known 
Instantons are self-dual classical solutions of the Yang- explicit instanton solutions on T 4 are quasi-abelian and 
Mills equations which play several important roles in of constant field strength 0. This lack of analytic in- 
Quantum Chromodynamics (QCD) [J H H H Q . Mucn formation gives additional motivation to the study of in- 
is known about their properties and their semiclassical stantons in lattice QCD, where instanton {i.e. self-dual) 
consequences in flat four-dimensional spacetime. On the configurations can be obtained by the technique of cool- 
other hand, at present the most powerful nonperturba- i n g, which has the goal of lowering the action S of the 
tive approach to QCD is lattice QCD, where spacetime configuration without changing its topological charge Q. 
is discretized. A discrete lattice with periodic boundary The action is bounded below by the topological charge 
conditions is in fact a four-toroidal mesh that approaches (in suitable units), and a numerical instanton configura- 
a continuum 4-torus in the continuum limit (a -» 0). tion is obtained when this bound is saturated. 
The properties of instantons on a continuum four-torus . 

T 4 differ in interesting ways from those of instantons on , The existence on T of instantons of topological charge 
R 4 @. In particular, while simple explicit expressions Ql > 2 was proved long ago by Taubes [10J. However, 

for periodic boundary conditions it can be proved that 
|Q| = 1 instantons cannot exist on T 4 This is an 

elegant corollary of a result known as the Nahm Trans- 
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ified in great detail [13( for the quasi-abelian constant 
field torus instantons. Nahm's construction has been 
generalized 0] to incorporate twisted boundary condi- 
tions, which describe additional non-abelian fluxes |15| . 
With such twisted boundary conditions it is possible to 
have a Q = 1 instanton solution, and powerful lattice 
techniques have been developed to numerically imple- 
ment the Nahm transform including twisted boundary 
conditions 0. This is a tour de force of lattice gauge 
theory, as it involves several highly nontrivial numeri- 
cal steps: cooling the configuration; finding zero modes 
of an associated Weyl-Dirac operator; reconstructing the 
Nahm transformed gauge field from the zero modes; and 
finally probing the self-duality properties of the trans- 
formed configuration. The results provide spectacular 
confirmation of the twisted Nahm transform construc- 
tion [H[l8|. 

In this paper we address a somewhat different ques- 
tion, relevant for instantons on an untwisted lattice. 
We accurately cool Monte-Carlo generated configura- 
tions in SU (3) with a highly improved action to produce 
highly self-dual configurations with a range of topological 
charges. Our goal is to probe in a detailed manner the 
difference(s) between the behavior under cooling of (un- 
twisted) configurations with topological charge \Q\ = 1 
and topological charge \Q\ = 2. The possibility of ob- 
serving the non-existence of self-dual \Q\ = 1 configu- 
rations in numerical simulations of lattice gauge theory 
on the 4-torus caught the attention of the lattice QCD 
community some time ago |19| . There, configurations 
with |Q| = I were prepared by cooling Monte-Carlo gen- 
erated ST/ (2) configurations with the standard Wilson- 
Gauge action [2(|. The \Q\ = 1 configurations were then 
cooled with improved actions toward the self-dual limit 
and were observed to be unstable to improved cooling. 
The topological structures were observed to shrink in size 
and eventually disappear from the lattice. This was pro- 
posed as lattice evidence for the non-existence of \Q\ = 1 
instantons on the four-torus [l9|. 

Here, we extend this analysis in several ways, with 
the goal of identifying and quantifying the impact of the 
Nahm transform corollary on \Q\ — 1 configurations in 
SU(3) gauge theory. We find three distinct signals of 
instability for the |Q| = I configurations. Interestingly 
these signals appear long before the would-be instantons 
have shrunk to the point where they "fall through the 
lattice". Furthermore, we show that these signals oc- 
cur only for the \Q\ — 1 configurations, and not for the 
individual separated lumps that make up our |Q| = 2 
configurations, which shows that this is really a global 
rather than a local effect, as expected for the continuum 
Nahm transform corollary [TT|. 

A key part of our analysis is a precise monitoring and 
minimizing of lattice discretization errors, which is cru- 
cial to interpreting the results of cooling studies. Er- 
rors in the discretization of the action can destabilize 
instantons without regard to the topological charge of 
the configuration. Improved actions are formulated to 



algebraically eliminate the leading finite lattice-spacing 
errors that arise as a result of approximating continuous 
space-time by a discrete mesh of points. This enables 
simulations to more accurately approach continuum be- 
havior while still retaining a finite lattice spacing |2*lj . 
The use of such improved actions in cooling algorithms 
has been shown to facilitate high-precision studies of the 
properties of lattice gauge fields [l£j, Hfl |22|, |23| ■ 

A convincing demonstration of the instability of \Q\ = 
1 lattice instantons requires careful monitoring of the in- 
stanton size and discretization errors. These are particu- 
larly important, as any instanton with a size smaller than 
the dislocation threshold of the cooling algorithm is elim- 
inated under improved cooling, again without regard to 
the topological charge of the configuration. The process 
is signified by large discretization errors associated with 
the small size of the object as it falls through the lattice. 

This point was recognized as a caveat in the interpreta- 
tion of the results of Ref. [n| , as the size of the topologi- 
cal structures in the configurations studied were initially 
very close to the dislocation threshold of the cooling al- 
gorithm. Thus, it was not possible to determine whether 
the disappearance of the structure was associated with 
the Nahm transform corollary, or simply the removal of a 
dislocation under improved cooling. Our results suggest 
that this disappearance stage of the topological charge 
evolution under cooling is in fact the final stage of an 
instability that manifests itself much earlier in the cool- 
ing process, when the topological objects are much larger 
than the dislocation threshold. 

Following from the discussion above, we propose the 
following criteria to identify convincingly the impact of 
the Nahm transform corollary: 

1. The lattice-discretized cooling action must be accu- 
rate with remaining errors acting to stabilize topo- 
logical structure. 

2. The action and topological charge density of a 
\Q\ = 1 configuration must be distributed over 
length scales much larger than the dislocation 
threshold of the improved cooling algorithm. 

3. Evidence must be presented to confirm that dis- 
cretization errors in the action arc minimal where 
the Nahm transform corollary manifests itself. 

4. Evidence that the action of the \Q\ = 1 configu- 
ration is approaching the self-dual limit of 8ir 2 /g 2 
must be provided. 

5. Moreover, one must demonstrate that the distribu- 
tion of action and topological charge distributions 
are approaching the classical form. 

We achieve each of these criteria, in the following manner: 

1. We utilize a highly- improved action free from Q (q 4 ) 
errors and with slightly positive 0(a e ) errors |22| 
to ensure the stability of instantons over several 
thousands of cooling sweeps. This action is known 



3 



as a three-loop improved action and appears better 
adapted to stabilizing topology than the five-loop 
improved action. 

2. The size of the action and topological charge dis- 
tributions are estimated by fitting the shape of the 
distribution surrounding the peak of the distribu- 
tion to the classical single-instanton density pro- 
file. The position and size, p, of the distributions 
are determined and compared with the dislocation 
threshold of the cooling action. 

3. Remaining lattice discretization errors are eval- 
uated by comparing the C(a 4 )-improved gauge- 
field action with an O(a 4 )-improved "reconstructed 
action" obtained from the square of an 0(a 4 )- 
improved lattice field-strength tensor 22], thus 
having different 0(a 6 ) discretization errors. 

4. Action and "reconstructed action" results will be 
reported in units of the single instanton action So = 

5. During the cooling process we also monitor the 
number of peaks in the action density. Even as 
the action approaches the single instanton action, 
numerous peaks in the action density can be iden- 
tified. Their position and size are also monitored 
as a function of cooling sweep. 

This paper is set out as follows. In Section [H] we 
briefly describe the highly-improved lattice discretiza- 
tion of the continuum action, field-strength tensor, re- 
constructed action, and topological charge operators. In 
Section IlIII we outline the simulation techniques and pa- 
rameters. Section llVl examines the evolution of \Q\ = 1 
configurations under improved cooling and puts the re- 
sults in perspective through a direct comparison with the 
behavior of \Q\ = 2 configurations. Our conclusions are 
presented in Section Ivl 

II. LATTICE ACTION AND TOPOLOGICAL 
CHARGE OPERATORS 

The lattice version of the Yang-Mills action was first 
proposed by Wilson [24[. The action is calculated from 

the plaquette, W$/ , a closed product of four link op- 
erators incorporating the link U^, 

Svm = PY,T,( 1 --fi Re trWji* 1 Hx)\ , (1) 

= \f d 4 xtrF^(x)+0(a 2 ), (2) 

provided (3 = 2N/g 2 for an SU(N) field theory. We will 

use the notation W^™ to denote the closed loop prod- 
uct (Wilson loop) in the \x — v plane with extent m lattice 
spacings in the /x-direction and n lattice spacings in the 



^-direction. Similarly, the lattice topological charge is 
obtained by summing the charge density over each lat- 
tice site, 

Q = 5Z^2 e ^ CT ^{Ff 1 v{x)F pa {x)} (3) 

x 

where p,, v, p, a sum over the directions of the lattice axes. 

Since the action and topological charge constructed 
from different Wilson loops have different 0(a 2 ) and 
higher terms we may cancel leading discretization errors 
by combining the contributions of different loops. For ex- 
ample, de Forcrand et al. 0,H(j have used tree-level im- 
provement to construct a lattice action which eliminates 
0(a 2 ) and 0(a A ) errors, by using combinations of up to 
five Wilson loop operators, which we denote Li....,L 5 . 
In the case of rectangular loops m ^ n the contribution 
of the loops in each direction is averaged, so that 

Li = W (lxl) 
I* = , 

L 3 = \{w^+W^}, 

L 5 = W#* 3 > . (4) 

The improved action of de Forcrand et al. is a linear com- 
bination of the Wilson actions calculated from L%, ...,£5, 
each divided by the relevant loop area squared (in units 
of a ), and respectively weighted by the constants 

ci = (19-55cb)/9, 

c 2 = (l-64c 5 )/9, 

c 3 = (640c 5 - 64)/45 , 

c 4 = 1/5 -2c5, (5) 

with C5 as a free variable, i.e., 

S = Cl 5(L0 + ^S(L 2 ) + ^S(L 3 ) + |s(£ 4 ) + ^S(L 5 ). 

(6) 

Selecting appropriate values for C5 enables one to set the 
contribution of certain loops to zero, creating so-called 
3-loop (c 5 = 1/10), or 4-loop (c 5 = 0) improved actions. 
Other values of C5 lead to 5-loop improved actions. Fol- 
lowing Ref. [l!J H3 , we consider c 5 = 1 /20 (midway be- 
tween the 3-loop and 4-loop values) to define the "stan- 
dard" 5-loop action. 

All of these actions have the same level of improve- 
ment, but in general they will have different 0(a e ) errors. 
We refer to the Wilson action, Eq. J5J), constructed from 
the plaquette as a 1-loop action. A 2-loop action may 
be constructed [2{| from |Li — \L 3 . These definitions 
of the variously-improved actions are used in our cooling 
algorithm. 

We construct a highly-improved lattice field-strength 
tensor, analogously to the construction of the cooling ac- 
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tion. The relations 



= {l + igjA.dx- 9 ^{j>A.dxf + 0{g^ , 

= {l-igj>A.dx~ 9 ^{j>A.dxf + 0{g 3 )Y 
indicate the field-strength tensor may be accessed via 

Y (V, xl) - ^ (lxl)t - 1* (w^ w£">t) 

= g <j> A.dx + 0{g 3 ) 



= ga 2 F^(x ) + 0(«?a 4 ) + 0(.gV) , 



(7) 



where we have subtracted one-third of the trace to en- 
force the traceless aspect of the Gell-Mann matrices. As 
with the improved action, the improved field-strength 
tensor is constructed from a combination of Wilson loops. 
Our improved field-strength tensor is given by 



(8) 



where Cj^' n ^ are the clover ave rag es ofmxn and n x 
m path-ordered link products |22(. The improvement 
constants take the values 





= 19/9- 


55fc 5 , 


k 2 


= 1/36- 


I6A5, 


k 3 


= 64fc 5 - 


32/45 , 


ki 


= 1/15- 


6fc 5 , 



and in this case k^ is the tunable free parameter. Fj^ p 
may be inserted directly into Eq. @ to create an im- 
proved topological charge operator, or into the equation 
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(9) 



to create a "reconstructed" improved action, Sr. Both 
improved operators will have errors of 0(a e ). Errors as- 
sociated with g are rapidly suppressed in the process of 
cooling. Appropriate choices of k 5 enable us to create 
3-loop or 4-loop improved versions of the field-strength. 
Our 5-loop improved field-strength tensor is defined with 
&5 = 1/180, midway between the 3-loop and 4-loop values 
of &5 = 1/90 and respectively While the 5-loop im- 
proved field-strength tensor is marginally more accurate, 
it is computationally expensive. We adopt the 3-loop 
improved field-strength tensor to construct the 0(a 4 )- 
improved topological charge and reconstructed action op- 
erators. 

Tadpole corrections to the improvement coefficients 
|25| are included and updated after every cooling sweep. 



Our experience |22j is that tadpole improvement factors 
arc beneficial in the early stages of cooling and remain 
beneficial even as uq — > 1. 

Previous investigations [2^] have indicated the 3-loop 
and 5-loop improved actions are similarly accurate in re- 
producing the classical action of approximately self-dual 
configurations. In some cases, the 5-loop action would 
underestimate the classical action. Such errors open 
the possibility of destabilizing topological structures over 
several thousand cooling sweeps, regardless of the topo- 
logical charge of the configuration. For the configurations 
investigated in Ref. |2!|] the errors in the 3-loop improved 
action were consistently positive; thus acting to stabilize 
topology. In the following, we adopt the 3-loop improved 
action. 



III. LATTICE APPROACH 

To investigate the stability of self-dual gauge field 
configurations on the (untwisted) 4-torus, we construct 
an ensemble of field configurations using the Cabibbo- 
Marinari [271 pseudo-heatbath algorithm with three diag- 
onal SU(2) subgroups looped over twice. We thermalizc 
for 5000 sweeps with an O(a 2 )-improved action from a 
cold start (all links set to the identit y) a nd select config- 
urations every 500 sweeps thereafter [2£| ■ Configurations 
are numbered consecutively in the order that they are 
produced in the Markov-chain process. Hence configura- 
tion 1 was saved after 5000 thermalization sweeps from 
a cold start, configuration 2 was saved 500 sweeps after 
configuration 1, and so on. Our results are generated 
on a 12 3 x 24 periodic lattice at (3 = 4.60, with a lattice 
spacing of a — 0.122(2) fm determined by a string tension 
analysis incorporating the lattice coulomb term. Cooling 
is performed with a 3-loop improved action with tadpole 
improvement (using the plaquette definition of the mean 
link) and the topological charge is assessed with the 3- 
loop improved operator. Links which may be updated in 
parallel are identified and partitioned [2!j . Updating the 
partitioned links in parallel minimizes the drift of objects 
in the cooled configurations. 

As a thermalized configuration is cooled, over the 
course of many hundreds of sweeps the action S of the 
configuration monotonically decreases. This occurs be- 
cause the cooling algorithm smooths out short-range fluc- 
tuations in the field. As the high-frequency components 
of the field are suppressed the underlying semi-classical 
structure of the field is revealed. If cooling proceeds for 
long enough the configurations will become self-dual con- 
sisting only of instantons or anti-instantons. 

The action, S, is bounded below by the magnitude of 
the topological charge, Q, i.e., 



S>S \Q\ = ^-\Q\ 

g 2 



(10) 



where in the continuum Q assumes an integer value 
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The quantity So = ^p- is the action associated with a 
single instanton and is independent of the instanton size. 

As cooling is applied uniformly over the 4-volume of 
the torus, configurations become locally self-dual. In 
the infinite volume limit where well-separated instantons 
and anti- instantons can simultaneously be present, the 
total action and topological charge will approach dis- 
crete values, satisfying the relations S/Sq = nj + tia 
and Q = m — tia, where n\ and ua are the number of 
instantons and anti-instantons present respectively. On 
a finite volume, one might expect to see some plateauing 
in the action density as a function of cooling sweep, when 
these conditions are satisfied. However, further cooling 
will lead to the annihilation of instanton-anti-instanton 
pairs in the finite volume, until only instantons or only 
anti-instantons remain. The configuration will then con- 
tinue to cool until it achieves complete global self-duality, 
i.e.,S/S = \Ql 

As discussed in the Introduction, self-dual configura- 
tions with \Q\ = 1 (i.e. a single instanton or a single 
anti-instanton) do not exist on the untwisted continuum 
four-torus Note, of course, that configurations with 
topological charge \Q\ = 1 do exist; they just cannot be 
made self-dual. In other words, the lower bound 5 > 5 
cannot be exactly saturated. We therefore expect to see 
different behavior of \Q\ — 1 and |Q| = 2 lattice config- 
urations under cooling. The next section is devoted to a 
study of signatures of this different behavior. 

IV. COOLING COMPARISON FOR Q = 1 AND 
Q = 2 CONFIGURATIONS 

A. Action and Charge Evolution 

Throughout the discussion of the results we will use the 
following notation: the cooling action will be denoted by 
5, the reconstructed action by Sr and the topological 
charge by Q. The type of improvement scheme used will 
be denoted by a number in parentheses. Hence the 2-loop 
improved cooling action, constructed from L\ and L% is 
written as 5(2), our topological charge calculated from a 
3- loop improved field-strength tensor is written as Q(3), 
our 3-loop reconstructed action (calculated from a 3-loop 
improved field-strength tensor) is written as Sr(3), and 
so on. The action reported throughout the cooling pro- 
cess will typically be divided by So to facilitate the com- 
parison with \Q\. Hence the curves in figures labeled as 
5(n) and Sr(ti) actually represent the normalized action 
values S(n)/So and Sn(n)/So respectively. 

1. \Q\ = 1 Configurations 

In Fig.^we show 5(3)/5q and |Q(3)| for configuration 
64. After a small number of cooling sweeps this config- 
uration rapidly settles down to a topological charge of 
one. Hence we expect that under continued cooling with 
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FIG. 1: Topological charge, Q(3), and normalized action, 
S(3)/So, for configuration 64. Notation is described in the 
text. In particular, Q and S are assessed here with the 3-loop 
improved topological charge and action operators. The action 
evolution flattens around 300 sweeps characteristic of a well 
separated instanton anti-instanton pair. Annihilation of the 
instanton anti-instanton pair follows as signified by the reduc- 
tion of the action by two units, while the topological charge 
remains unchanged. The nearly self-dual single instanton con- 
figuration destabilizes around sweep 2700. 



negligible discretization errors, it should never reach true 
self-duality. We do observe a long plateau where the con- 
figuration appears to closely approach self-duality. How- 
ever, around cooling sweep 2700 it eventually destabilizes 
and collapses to triviality. This behavior is consistent 
with earlier results 0] for SU(2). We will argue in the 
following that this final destabilization of the instanton is 
in fact the final stage of a series of evolution stages, with 
this final disappearance of the would-be instanton being 
due to its small size relative to the dislocation threshold 
of the cooling algorithm. 

Let us now examine the same configuration with an 
enhanced vertical scale as in Fig. |2 The configuration 
fails to achieve self-duality before it destabilizes and be- 
comes trivial. However, a very high level of accuracy is 
obtained for the topological charge, Q, between sweeps 
200-400. This is the region in which the action evolution 
flattens with S/Sq ~ 3 signaling approximate local self- 
duality. This region provides a benchmark for the level of 
accuracy that may be achieved through the use of highly 
improved operators. 

For sweeps 400-1400 we see the value of Q rise away 
from the integer value of one. This "cresting behavior" 
commences at the same time that the action leaves the 
regime of S/Sq ~ 3 and approaches the classical instan- 
ton action. We note that it is well separated from the 
ultimate collapse of the configuration. 

The two 3-loop improved actions (5(3) and 5^(3)) re- 
main greater than 5o|Q| as required. The two actions 
deviate significantly from one another when the configu- 
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ration collapses to triviality, which is indicative of large 
discretization errors. Thus, this destabilization, at this 
stage of the cooling process, should not be attributed to 
the Nahm transform corollary. In contrast, the two ac- 
tion measures agree relatively well even at the peak in 
the topological charge cresting. Note that S(3)/Sq > 1, 
so that errors should act to stabilize the configuration. 

To illustrate the reproducible nature of these results we 
also present results in Fig.|31for configuration 11; another 
|Q| = 1 configuration. Beyond 200 sweeps, the behavior 
is nearly identical. Beyond approximately 1500 cooling 
sweeps we see a crossing followed by a large divergence 
of the action values and a sudden drop in the topological 
charge. This indicates the onset of severe discretization 
errors, again suggesting a small size for the nontrivial 
topological structure, the "would-be" instanton. 

Based on these results, we propose that the cresting of 
the topological charge early in the cooling process, and 
associated with the departure of 5(3) / Sq from the regime 
of 3 toward 1, is a distinct signature of the ultimate insta- 
bility of the single instanton on the lattice. The eventual 
disappearance of the would-be instanton happens much 
later in the cooling process, and appears to rather be a 
consequence of large discretization errors in that region. 
In the next sub-section, we contrast this directly with the 
behavior of \Q\ = 2 configurations under cooling. 



2. \Q\ = 2 Configurations 

To show that this cresting behavior is specific to \Q\ = 
1 configurations, we compare one and two instanton con- 
figurations in Figs. 01 and [5] In these figures we consider 
the three-loop quantities 5(3)/5o and |Q(3)|. For the 
two-instanton case, we plot (5(3)/5o) — 1 and |Q(3)| — 1, 
so that the curves for the \Q\ = 1 and \Q\ = 2 cases may 
be directly compared. The level of self-duality achieved 
by the two-instanton configurations is far superior to that 
of the single-instanton configurations. The cresting be- 
havior observed in the topological charge of the \Q\ = 1 
configurations is not seen at all for the \Q\ = 2 configu- 
rations. 

We stress that the behavior of the \Q\ = 1 configura- 
tions differ significantly from that of the \Q\ = 2 configu- 
rations long before the \Q\ — 1 configurations collapse to 
triviality. Indeed, this deviation becomes manifest rela- 
tively early in the cooling process. For configuration 64 
of Figs . ^ and [21 the cresting behavior begins as the num- 
ber of locally self-dual objects on the lattice drops from 
three to one. 

To confirm the consistency of this picture, it is im- 
portant to know the size of these objects in order to 
evaluate whether the dislocation threshold of the cool- 
ing algorithm is playing a role in the evolution of the 
configurations, and at what stage of the cooling process 
this is important. This is the subject of the next section. 
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FIG. 2: Plot showing the fine detail of the \Q\ = 1 configu- 
ration illustrated in Fig. configuration 64. The three-loop 
cooling action, 5(3), three- loop reconstructed action, Sr(3), 
and the three-loop topological charge, Q(3) are plotted. A 
very high level of accuracy is obtained for the topological 
charge between sweeps 200-400. A "cresting behavior" ap- 
pears as S/So drops from ~ 3 toward 1. 
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FIG. 3: The fine detail of another \Q\ — 1 configuration; 
configuration 11. Beyond 200 sweeps, the behavior is nearly 
identical to that illustrated in Fig. [5] 



B. Instanton size evolution 

To investigate further the behavior of the topological 
structure of these |Q| = 1 and \Q\ = 2 configurations 
under cooling, we implement an algorithm which iden- 
tifies local peaks in the action and topological charge 
densities. Peaks in the action density are identified by 
finding a point at the center of a 3 4 hypercube whose 
action density exceeds that of the neighboring 80 points 
of the hypercube. The algorithm may also be applied to 
the topological charge density in two steps, reversing the 
sign of the topological charge density to convert valleys 
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FIG. 4: Comparison of the \Q\ = 1 configuration (configu- 
ration 64, solid line) with a \Q\ — 2 configuration (configu- 
ration 90, dashed line). For the two-instanton configuration 
both S(3)/So and \Q\ have been reduced by unity so that the 
curves may be directly compared with the \Q\ = 1 configura- 
tion. The accurate behavior of the two-instanton configura- 
tion contrasts the single-instanton case. 
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FIG. 5: A similar comparison of a \Q\ = 1 configuration (con- 
figuration 11, solid line) with a \Q\ = 2 configuration (config- 
uration 27, dashed line). As in Fig. S(3)/S and \Q\ for 
the two-instanton configuration have been reduced by unity 
to allow a direct comparison of the curves with the \Q\ = 1 
configuration. 



to peaks. 

The algorithm then fits the structure of the gauge fields 
around these peaks to the classical instanton form for the 
action density 



S(x) = ^jj— 

7T Z ((X 



xq) 



p 2 )Y 



(11) 



generalized by the inclusion of £ to allow for an over- 
all normalization different from 1 due to periodic images 



of the action density |30j. For an instanton in infinite 
volume on R 4 , £ = 1. The fit parameters returned by 
this algorithm are the coordinates of the center of each 
instanton- like peak, xq, the instanton-size parameter for 
each peak, p, and the overall scale factor for the peak, £. 

It should be noted that not every peak found by this 
algorithm is an instanton, particularly during the early 
stages of the cooling process. However, we certainly an- 
ticipate that the peaks which survive under cooling will 
approach the form of instantons or anti-instantons as 
the cooling proceeds. Peaks which do not correspond to 
instanton-like structures will disappear in the early stages 
of cooling. Of the larger would-be instantons and anti- 
instantons surviving these early cooling steps, all but \Q\ 
of them will ultimately disappear due to the annihilation 
of instanton-anti- instanton pairs. 

We use this algorithm to investigate the topological 
properties of configurations 1 1 and 64 (corresponding to 
\Q\ = 1) and in configurations 27 and 90 (corresponding 
to \Q\ — 2). In Fig. we illustrate how the instanton 
size, p, varies with sweep number for configurations 64 
and 11. Three initial instanton-like peaks are identified 
in configuration 64 and four are identified in configura- 
tion 11. The solid line in each of these figures denotes 
the peak which survives long-term cooling and becomes 
the approximately self-dual object in each configuration. 
The dashed lines denote the temporary peaks, which dis- 
appear under cooling. 

In configuration 64 we see that all peaks rapidly grow 
in size initially. Both of the temporary peaks are seen to 
expand in size with the associated peak height decreasing 
until the peaks "melt" away. These peaks are associated 
with the plateau in the action observed in Fig. ^ Hence 
we see that topological objects can disappear on the lat- 
tice by "melting away" where p — > oo, which is a different 
mechanism from the behavior of "falling through the lat- 
tice" where p — > 0. 

The remaining topological object reacts by slowing its 
growth. Once all of the other structures have disap- 
peared, the surviving object shrinks rapidly. In minimiz- 
ing its size relative to the finite volume of the 4-torus, 
the object can best approximate the self-dual nature of 
an instanton and temporarily evade the ultimate conse- 
quence of the Nahm transform. The rate of shrinkage is 
related to the size of the object. The shrinkage slowly 
continues until the instanton size is approximately two 
lattice spacings, at which time the object is smaller than 
the dislocation threshold of the cooling algorithm and is 
rapidly suppressed. 

Similar results are observed for configuration 11 (an- 
other \Q\ = 1 configuration), as shown in the lower plot of 
Fig. Again all peaks rapidly grow in size initially, but 
as each of the three temporary peaks expands and "melts 
away" , the single instanton reacts by slowing its growth. 
Once all of the other structures have disappeared, we see 
that the surviving single (anti-)instanton shrinks rapidly 
due to its large size. The shrinkage continues until the 
dislocation threshold is encountered, at which point it is 
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FIG. 6: Size of the single instanton (p measured in lattice 
units) in configuration 64 (upper figure) and configuration 
11 (lower figure). An abrupt change in the size evolution 
occurs as other peaks in the distribution disappear through 
the process of "melting" as described in the text. After sweep 
2500 the instanton in each configuration rapidly disappears 
and hence the associated size parameter cannot be measured 
accurately. Dashed lines represent the size of transitory peaks. 



rapidly suppressed. 

The point of inflection seen in the evolution of the 
single-instanton size around 1900 sweeps in Fig. signi- 
fies the dislocation threshold of the S(3) improved cool- 
ing algorithm. For both configurations p ~ 2.2 at the 
point of inflection. This dislocation threshold is similar 
to po = 2.23 for S(5) as determined in Ref. |19j . 

In the continuum infinite volume limit, instantons have 
no implicit scale and can be any size with the same ac- 
tion Sq. On a finite 4-torus in the continuum, the only 
relevant quantity is the relative size of the instanton to 
the 4-torus size. On the continuum 4-torus, we therefore 
expect that in the limit where the size of the would-be in- 
stanton vanishes with respect to the 4-torus size, it should 
able to approach arbitrarily closely to self-duality. Thus, 



under continued cooling on the lattice a \Q\ = 1 con- 
figuration will approach the self-dual limit by decreasing 
its size. As it does so, discretization errors will become 
increasingly large until the object encounters the disloca- 
tion threshold, at which point the object will be removed 
by the cooling algorithm. This is the ultimate fate of a 
|Q| = 1 configuration under cooling. 

In both \Q\ = 1 cases, the onset of the cresting in 
the topological charge occurs when the size of the object 
is seven to eight lattice spacings, well above the size of 
two lattice spacings where the dislocation threshold of 
the cooling algorithm removes the objects. Thus, at the 
onset of the cresting of \Q\ the objects are large enough 
that the dislocation threshold of the cooling algorithm is 
irrelevant. Since the cresting behavior appears to be a 
distinct signal of the ultimate disappearance of the |Q| = 
1 objects, this suggests that the dislocation threshold is 
not the primary reason for the instability of the \Q\ = 1 
configurations, even though it is the reason for its final 
disappearance. 

To make this even more concrete, we contrast these re- 
sults for the size dependence of the |Q| = 1 objects under 
cooling, with the corresponding results for the \Q\ = 2 
configurations. 

The sizes of the two instantons in each of configura- 
tions 27 and 90 (which each have |Q| = 2) are shown in 
Fig. [7| First, note that in each case the configuration 
consists of two identifiably separate topological objects, 
rather than one single object of instanton number 2. For 
each configuration, the two instantons drift apart under 
cooling, appearing to repel one another. More interest- 
ingly, in each case, the size of one instanton grows while 
the other shrinks, until they each settle (after roughly 
1000 sweeps) to approximately the same size. We do not 
know if this behavior is necessary - however, the fact 
that the same thing happened in these two independent 
cases is suggestive. Perhaps this symmetry between the 
sizes of the cooled two-instanton configuration is due to 
an image effect that balances two objects at equal sizes 
roughly at "opposite ends" of the lattice. It would be in- 
teresting to probe this question more deeply, both on the 
lattice and for the continuum four-torus. For example, 
for SU(2) ADHM instantons on R 4 it is known that the 
behavior of Q = 2 instantons depends crucially on the 
relative SU(2) orientation of the two constituent instan- 
tons, which has important implications for the instanton 
size distribution |31| . 

Another striking contrast between the size-dependence 
of the |Q| = 2 and \Q\ — 1 configurations is that the 
\Q\ — 1 objects shrink in size by a factor of approxi- 
mately three as the cooling proceeds over the course of 
two thousand sweeps, while in the \Q\ = 2 configurations, 
each object changes size by no more than 10%. This in- 
dicates a significant level of stability throughout most 
of the course of the cooling process. Since the |Q| = 2 
configurations consist of two distinct topological objects, 
each of which may be thought of as a single instanton, 
it is striking that for these configurations the single in- 
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FIG. 7: Size evolution of the two instantons (p measured in 
lattice units) in configuration 27, and the two instantons in 
configuration 90. The solid and dashed lines represent the 
size of one of the two instantons present in each configura- 
tion. The discontinuity in the lower figure (configuration 27) 
is believed to be caused by a sudden change in parameter 
values as the fitting program used to calculate the single in- 
stanton size from our lattice data jumped from one local x 2 
minimum to another. 



stantons do not shrink in the same way that the \Q\ = 1 
would-be instantons do. This shows that the distinction 
between the \Q\ = 1 and \Q\ — 2 configurations is really 
a global effect, not a local one, and this is exactly what 
we would expect from the Nahm transform corollary |llj| . 
So, this difference in the size dependence under cooling 
is another signal (in addition to the cresting behavior of 
\Q\ identified in the previous section) of the instability of 
|Q| = 1 configurations under cooling, and once again we 
note that the difference can be seen already very early 
in the cooling process, rather than just at the late time 
when the \Q\ = 1 objects vanish due to discretization 
errors. 



C. Discretization Errors 

In this last section we identify a third distinct signal of 
the instability of the \Q\ — 1 configurations under cool- 
ing, and show that it can also be seen early in the cooling 
process. It was noted already in Figs. [21 and [21 that S and 
Sr diverge rapidly in the period prior to the collapse of 
the \Q\ — 1 configurations. Actually, for the \Q\ = 1 
configurations, the deviation between these two actions 
appears much earlier, as is shown in Fig. [SI In Fig. [SJwc 
present plots of the ratio of S to Sr, for both \Q\ = 1 
configurations (11 and 64) and \Q\ = 2 configurations (27 
and 90). Since the cooling action and reconstructed ac- 
tion are improved differently, they are expected to have 
different 0(a 6 ) discretization errors. The ratio of S to Sr 
therefore indicates the relative scale of the 0(a 6 ) errors 
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FIG. 8: Ratio of the cooling action to the reconstructed action 
for configurations 64 and 90 (top) and 11 and 27 (bottom). 
Configurations 11 and 64 are \Q\ ~ 1 configurations, whereas 
configurations 27 and 90 are \Q\ = 2 configurations. 



present. 

Early in the cooling evolution, peaks in S/ Sr are seen 
to be associated with the disappearance with topolog- 
ical structures as illustrated in Fig. In the case of 
the single-instanton configurations, the ratio of the cool- 
ing action to reconstructed action, S/Sr, diverges from 
unity slowly during the middle stages of the cooling pro- 
cess. The point where this ratio starts falling below 1 
corresponds to the beginning of the cresting behavior of 
\Q\ which was identified for the \Q\ = 1 configurations in 
Section llV Al . This corresponds to the early stage where 
the difference between the size evolutions of the |Q| = 1 
and \Q\ = 2 configurations becomes clear. Also, the po- 
sition of the valley in the ratio S/Sr corresponds to the 
peak of the cresting of the topological charge illustrated 
in Figs. H and 

Finally, S/Sr diverges rapidly as the object present 
in each configuration falls through the lattice. As the 
instanton shrinks, the 0(a 6 ) errors overwhelm the im- 
proved action operator and large errors allow the disap- 
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pearance of the object. By contrast, in the two \Q\ = 2 
configurations the ratio S/Sr remains within 0.02% of 
unity indefinitely. 



V. CONCLUSIONS 

In this paper we have compared the cooling of \Q\ = 1 
and \Q\ = 2 configurations in SU(3) gauge theory, us- 
ing improved lattice actions and precise monitoring of 
the sizes of the topological objects and of the lattice dis- 
cretization errors. In the continuum it is not possible to 
have a self-dual \Q\ = 1 configuration on the untwisted 
four-torus [11,], and the results from our lattice analysis 
suggests three distinct signals, each of which occurs very 
early in the cooling process, that the \Q\ = 1 configu- 
rations behave differently from the \Q\ = 2 ones, and 
that they will ultimately shrink to below the discretiza- 
tion threshold and then disappear much later in the cool- 
ing process. These instability signals occur at a stage in 
the cooling when the topological objects are still much 
larger than the lattice discretization threshold. None of 
these signals is present for the \Q\ = 2 configurations we 
studied, even though these consisted of two isolated sin- 
gle instantons. This fact is a clear indication that these 
|Q| = 1 instability signals are reflecting the global prop- 
erties of the torus, which is the essence of the continuum 
torus result [ll|, rather than local effects such as those 
that ultimately cause a very small |Q| = 1 object to fall 
through the lattice due to discretization errors. 

The three instability signals we have found to be char- 
acteristic of the \Q\ = 1 configurations are as follows. 
First, there is a cresting behavior of the topological 
charge \Q\ away from its precise integer value of 1, as 
shown in Figures 2 and 3. This does not happen for 
the \Q\ = 2 configurations, as shown in Figures 4 and 5. 
Second, for the |Q| = 1 configurations there is an initial 
swelling of the size of the would-be instanton, but then 
a steady shrinkage begins at the same point where the 
cresting in |Q| is observed. This is illustrated in Figure 
6. Once again, this is completely different from the be- 
havior for the \Q\ = 2 configurations, which is shown in 
Figure 7. The third signal is that for the \Q\ = 1 con- 
figurations the ratio S/Sr of the cooling action to the 
reconstructed action begins a steady fall below one, be- 
fore bottoming out and eventually diverging much later 
in the cooling process, as shown in Figure 8. This devi- 
ation below a ratio of 1 begins at the same point, early 
in the cooling process, where the cresting of \Q\ and the 
rapid shrinkage of the object begin. This is also very dif- 
ferent from the behavior of the \Q\ = 2 configurations, 
already at this early stage where the topological objects 
are large. 

These results suggest the following picture of the insta- 
bility, and ultimate disappearance, of \Q\ = 1 configura- 
tions under cooling. The first stage sees the rapid growth 
and disappearance of temporary peak structures. The 
second stage, whose onset is indicated by the three sig- 
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FIG. 9: The behavior of instanton peak size, p (upper solid 
line, right scale), against topological charge (lower solid line, 
left scale) and Sr — S + l (dashed line, left scale), as a function 
of sweep number, for \Q\ ~ 1 configurations 64 (upper figure) 
and 11 (lower figure). Notice that Sr — 5 + 1 and Q move 
away from integer values at the same time that the instanton 
starts to shrink. Cooling is performed with the 5(3) action, 
and the topological charge is Q(3). We have taken the value 
Sr — 5 + 1 rather than Sr — S so that the action and charge 
curves can be plotted on the same precise scale. 



nals listed in the previous paragraph, sees the shrinkage 
with cooling of the single instanton. This type of shrink- 
age does not happen for the separate constituents of the 
\Q\ = 2 configurations. During this shrinkage stage the 
topological charge deviates slightly from its integer value 
of 1. The third and final stage is the rapid disappear- 
ance of the single instantons due to dislocation errors at 
an instanton size of approximately two lattice spacings. 

To emphasize these points, we show in Fig. [21 the in- 
stanton size p, the topological charge Q, and the differ- 
ence Sr — 5+1 between the reconstructed action and 
cooling action for the | Q | = 1 configurations 64 and 1 1 in 
the region where the Nahm transform corollary manifests 
itself. (The addition of one unit to the difference serves 



11 



to plot the action and charge curves on the same precise 
scale.) Since «Sr and S have different (D(a 6 ) errors, their 
difference is an extremely effective probe of the scale of 
discretization errors in the configuration. As the instan- 
ton changes from expanding behavior to shrinking behav- 
ior, the topological charge and difference Sr — S + 1 are 
both almost identically equal to one, indicating minimal 
discretization errors at the onset of the Nahm transform 
corollary. This strongly suggests that the sudden change 
in behavior of the would-be instanton is due to the global 
nature of the field on the toroidal lattice, and not to local 
discretization errors. 

The |Q| = 2 configurations behave very differently 
from the \Q\ = 1 configurations. For each of the two 
|Q| = 2 configurations we studied, the configuration con- 
sisted of two isolated lumps that drift apart under cooling 
and have the remarkable property that one lump shrinks, 
while the other swells, until they reach the same size. It 
would be interesting to know if this behavior is generic, 
due to some torus periodicity-induced balance, and if 
there is any trace of such behavior in the continuum. Un- 
fortunately, this final question is a difficult one, as there 
are no known nontrivial (i.e. inhomogeneous) analytic 
torus instantons. 

We close by noting that all five criteria outlined in the 
introduction have been met. In particular: 

1. The S(3) cooling action is very accurate with re- 
maining errors positive, acting to stabilize topolog- 
ical structure. 

2. The action and topological charge density of the 
|Q| ~ 1 configurations are seen to be distributed 



over large length scales of 7 to 8 lattice spacings, 
much larger than the dislocation threshold of the 
improved cooling algorithm at 2.2 lattice spacings. 

3. Comparisons of the reconstructed action with the 
more traditionally constructed improved action 
have provided a powerful way of investigating the 
scale of discretization errors in the configurations 
under investigation. 

4. |Q| = 1 configurations have been shown to be stable 
until the action S drops below 3 So and approaches 
the self-dual limit of Sq = 8ir 2 / ' g 1 . Extra peaks in 
the action density disappear as S/Sq drops from 3 
toward 1. 

5. Finally, the action and topological charge distribu- 
tions are seen to approach the classical instanton 
form with only a single peak appearing and the in- 
stanton scale parameter £ — > 1 as cooling proceeds. 
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